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SIMPLODE: 

AN  IMPLODING  GAS  PUFF  PLASMA  MODEL 
L  NEON 


I.  INTRODUCTION 

The  use  of  a  gas  puff  plasma  as  a  transducer  for  the 
conversion  of  electrical  to  radiative  energy  has  been  the  focus 
of  an  intense  experimental  and  theoretical  program  for  the  past 
several  years.  Although  a  considerable  amount  of  electrical  and 
spectral  data  has  been  generated  the  bulk  of  it  remains 
uncorrelated  and  unanalyzed.  Only  for  those  shots  where  the 
radiative  yield  was  considered  impressive  has  there  been  any 
analysis.  Generally,  the  "analysis"  provided  information  on  the 
average  plasma  temperature  and  density  from  emission  spectra, 
plasma  size  from  x-ray  pinhole  photos  and  total  radiative  yield 
from  calorimetry.  Based  on  this  "analysis"  it  would  be  difficult 
to  obtain  anything  but  the  most  rudimentary  understanding  of  the 
plasma  dynamics,  certainly  not  enough  to  allow  the  systematic 
selection  of  both  the  machine  parameters  and  load  characteristics 
that  would  predict  an  optimum  radiative  yield  shot. 

The  theory  effort  is  chartered  to  develop  a  variety  of 
theoretical  models,  tools,  and  numerical  simulations  to 
understand,  guide,  and  predict  the  behavior  of  these  high 
brightness  laboratory  x-ray  sources.  The  models  vary  from  the 
simple  0-D  hydrodynamic  implosion  schemes  to  the  complex  and 
sophisticated  2-D  Non-LTE  radiation  magnetohydrodynami c  numerical 
simulations.  These  models  have  recently  been  used  successfully 
for  analysis  and  interpretation  of  several  past  and  present 
experiments  as  well  as  for  scaling  to  higher  energy  systems.1-^ 
For  this  investigation  it  was  essential  to  have  a  theoretical 
capability  that  would  provide,  in  real  time,  an  interactive  link 
between  theory  and  experiment.  This  requirement  precluded  the 
U3e  of  the  more  sophisticated  numerical  models  which  required 
large  continuous  intervals  of  computer  time,  which  was 
unavailable,  and  therefore  dictated  the  application  of  a  less 
Manuscript  approved  May  1, 1986. 


sophisticated  model  but  one  that  included  some  of  the  relevant 
physics.  The  SIMPLODE  model  satisfies  these  requirements. 
Briefly,  the  SIMPLODE  model  describes  the  0-D  implosion  dynamics 
of  a  current  driven  plasma  slug  self-consistently  coupled  to  a 
Non-LTE  radiation  physics  model  and  is  ideal  for  use  with  the  gas 
puff  experiments  at  NRL  on  the  GAMBLE  II  facility. 

Recently  the  GAMBLE  II  pulse  power  facility  has  been 
upgraded  to  accomodate  gas  puff  loads.  This  modification 
enhances  GAMBLE  II’s  versatility  by  expanding  the  types  of 
material  loads  that  can  be  investigated  on  it  making  the  facility 
very  attractive  for  examining  Plasma  Radiation  Source  (PRS) 
Loads.  The  inaugural  series  of  experiments  involved  krypton, 
argon,  and  neon  with  most  of  the  work  being  devoted  to  neon. 
These  experiments  were  designed  to  calibrate  the  facility  and 
determine  its  performance  level  in  dealing  with  gas  puff  loads 
and  nozzle  technology.  Also,  neon  was  a  better  matched  load  for 
GAMBLE  II  in  terms  of  the  efficient  conversion  of  electrical 
energy  to  K-shell  radiation.  Another  aspect  of  the  GAMBLE  II 
experiments  was  to  examine  the  impact  of  sharper  current  risetime 
on  the  implosion  dynamics  and  radiative  yields. 

The  work  presented  here  will  be  concerned  primarily  with 
some  of  the  theoretical  aspects  of  the  neon  gas  puff  plasmas.  In 
particular,  the  dynamic  and  radiative  characteristics  of  the 
implosion  as  a  function  of  input  parameters.  A  comparison  with 
some  of  the  experimental  results  will  be  made.  A  report 
describing  the  majority  of  the  experiments  is  currently  in 
preparation  in  conjunction  with  the  Plasma  Technology  Branch  and 
is  forthcoming. 

I I .  Model 

The  simplest  description  of  an  imploding  Z-pinch  plasma  is 
probably  the  Bennet  pinch  equilibrium  model.1*  It  represents  an 
equilibrium  balance  between  the  fluid  and  the  magnetic  pressures 
in  combination  with  an  equilibrium  balance  between  the  energy  in 
and  the  energy  out.  The  obvious  extension  to  the  Bennet 


equilibrium  pinch  model  is  the  inclusion  of  flow  parameters 
following  the  temporal  evolution  of  the  plasma,  i.e.,  maintain 
the  simple  philosophy  of  the  Bennet  pinch  but  allow  the  plasma  to 
evolve  in  time.  In  essence,  this  is  our  philosophy  -  a  dynamic 
Bennet  pinch.  Also,  like  Shearer,^  who  included  a  radiation 
cooling  term  in  the  form  of  bremsstr ahl ung  losses,  we  also 
include  radiation  cooling  but  in  a  more  extensive  fashion.  The 
radiation  cooling  term  in  our  model  includes  contributions  from 
free-free,  free-bound,  and  bound-bound  transitions  and  are 


determined  from  a 
ionization  dynamics. 


Non-LTE  co 1 1 i s i onal -r ad i at i ve  model  of  the 
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The  model  -  SIMPLODE  -  treats  a  cylindrical  annular  gas  puff 
plasma  of  uniform  density  carrying  a  uniform  current  in  the  Z- 
direction.  Only  the  radial  motion  is  considered  -  hence,  the 
plasma  is  always  uniform  in  the  Z-direction,  i.e.,  no  axial 
structure.  The  radial  acceleration  is  determined  from  the  force 
equation ,  vig . 
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where  ly^  is  essentially  the  sum  of  ionization  energies  and  is 

sometimes  loosely  referred  to  as  chemical  potential.  Acurr  “ 
2  2 

it  ( R  -  R,  )  where  R.(R<)  represent  the  outer  (inner)  radius 

°  .  .  ra/rai 

of  the  plasma,  V  «  2  i  l  r  r,  N  =  — where  m  is  the  mass,  m^ 

irr  l 

the  atomic  mass  and  the  thermal  pressure  is  N(1+z)kT.  The  first 
term  on  the  RHS  of  Eq.  (2)  represents  the  work  done  in  com¬ 
pressing  the  plasma,  the  second  term  is  the  joule  heating  source, 
and  the  third  term  is  the  power  radiated  per  unit  volume  from  the 
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volume  V.  n  represents  classical  resistivity,  and  l  the  plasma 
length.  The  self-consistent  solution  to  these  equations  is  what 
we  refer  to  as  SIMPLODE.  In  the  absence  of  a  circuit  model,  I  ( « 
current  waveform)  is  determined  experimentally  and  used  as  an 
input  parameter.  [The  model  is  presently  being  extended  to 
evaluate  s  el  f- co  n  3  i  s  t  entl  y  the  current  from  a  circuit  model  based 
on  the  open  circuit  voltage.]  It  i3  not  yet  clear  that  there  are 
any  significant  differences  between  the  current  waveform  as 
recorded  from  the  current  monitors  and  the  theoretically 
calculated  current  waveform,  except  at  late  time,  i  .e  .  ,  well 
after  peak  current,  where  the  current  crowbars  and  does  not 
vanish  on  the  experimental  waveform.  If  taken  seriously,  it  can 
maintain  the  plasma  in  a  heated  compressed  state  and  overestimate 
the  radiative  losses,  i  .e  .  ,  predict  a  higher  yield.  Hence,  the 

input  current  waveforms  have  been  reduced  to  zero  at  late  time  to 
avoid  over  es  timating  the  yield. 

III.  Results  and  Discussion 

Since  most  of  the  experiments  were  done  with  neon  ga3  puffs, 
the  numerical  simulations  will  also  focus  on  neon.  In  addition, 
the  experiments  were  divided  into  two  categories:  with  and 

without  the  P_lasma  Erosion  Opening  Switch  (PEOS).  To  differ¬ 
entiate  between  the  categories  we  will  generally  refer  to  them 
as  the  switch  and  no-switch  cases,  respectively.  However,  when 
the  PEOS  was  included  both  the  peak  current  and  current  waveform 
were  substantially  different  from  the  no-’switch  case.  In  fact, 
the  peak  current  only  achieved  70%  of  the  no-switch  value  and 

exhibited  a  much  broader  current  waveform.  These  differences 
will  be  quantitatively  delineated  below. 

The  SIMPLODE  code  requires  as  input  data  the  current 
waveform,  mass  per  unit  length,  and  the  inner  and  outer  radius  of 
the  plasma  slug.  The  plasma  length  is  taken  to  be  4  cm  in  all 

the  calculations  presented  here.  For  different  lengths  the  total 
yield  can  be  appropriately  scaled.  The  selected  cases  considered 


are  presented  in  Table  1. 


TABLE  1 


j  ( ugm/ cm  ) 

RB ( cm  ) 

R i (cm  ) 

23.5 

1  .45 

1  .  05 

23.5 

1  .  55 

0.55 

No-Swi tch 

23.5 

1  .  786 

0.655 

25.0 

1  .  55 

0.95 

35.0 

1  .  55 

0.95 

.23.5 

1  .45 

1  .  05 

Switch 

23-5 

1  .  55 

0.95 

23-5 

1  .  786 

0.655 

35.0 

1  .  55 

0.95 

The  experimental  current  waveform  for  the  no-switch  cases  is 
shown  in  Fig.  1.  The  current  reaches  a  peak  value  of  1.25  MA, 

has  a  current  ri3etime  of  about  40  ns,  and  a  FWHM  of  90  ns .  The 

measured  current  waveform  as  recorded  by  the  current  monitors 
crowbars  at  about  1  30  ns  and  maintains  a  steady  constant  value. 
In  reality,  the  current  vanishes  at  late  time.  We  have  taken  the 
liberty  of  extending  it  to  zero  at  135  ns  to  insure  that  it  does 
vanish.  This  current  waveform  has  been  used  for  all  the  no¬ 

switch  simulations.  This  is  strictly  not  correct  since  the  load 
impedance  varies  in  time  and  this  is  not  reflected  in  a  modified 
current  trace.  In  the  future,  SIMPLODE  will  be  coupled  to  a 
circuit  model  that  uses  the  open  circuit  voltage  to  self- 
consistently  predict  a  current  for  a  time  varying  load.  As  an 
example  of  the  no-3witch  case  we  present  the  results  of  the 
SIMPLODE  code  for  the  current  trace  shown  in  Fig.  1  and  the 

initial  conditions:  M/4  »  30  pgm/cm,  RQ»1.55  cm  and  Rj-0.95  cm. 
Figure  2  shows  the  time  evolution  of  the  inner  and  outer 
radius.  Note  that  the  outer  radius  pinches  down  to  about  0.5  mm 
and  bounces  back  out  at  150  ns,  i.e.,  the  fluid  pressure  exceeds 
the  magnetic  pressure  and  pushes  the  plasma  out.  This  happens 
s  i  gn  i  f i cantl y  after  peak  current.  The  implosion  velocity  of  the 
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outer  surface  is  shown  in  Fig.  3  and  reaches  a  maximum  value  of 
2.1  cm/ysec  and  is  zero  at  the  bounce  time.  The  effective  charge 
state,  ZEFF,  is  shown  in  Fig.  4  and  is  merely  a  reflection  of  the 
temperature  and  density  of  the  plasma.  The  sharp  spikes  are  not 
real  but  are  numerical  artifacts  associated  with  the  time-step 
used  in  this  particular  calculation.  If  a  significantly  smaller 
time-step  was  used  these  artifacts  would  disappear  leaving  a 
smooth  curve.  However,  reduction  in  At  leads  to  significant 
increases  in  computer  time.  Test  runs  with  smaller  At’s  produce 
values  equivalent  to  "eyeball"  averaging  and  smoothing  of  the 
results  presented  in  Fig.  4.  The  effective  charge  state  of  the 
plasma  is  about  3  at  50  ns.  The  temperature  is  shown  in  Fig.  5 
and  reaches  a  peak  value  of  about  125  eV  around  140  ns.  This 
occurs  well  after  peak  current  and  during  the  final  compression 
phase  before  bouncing  out.  The  ion  density  peaks  around  150  ns 
with  a  value  of  2.2  X  102®  cm~3  as  shown  in  Fig.  6.  Both 
temperature  and  density  peak  during  the  final  compress ional 
phase.  The  ion  density  has  a  narrow,  well  peaked  delta-function 
like  profile  representative  of  a  good  compressional  implosion. 
The  local  K-shell  radiative  yield  is  shown  as  a  function  of  time 
in  Fig.  7.  A  sharp  pulse  is  emitted  over  12  ns  beginning  at 
about  135  ns  and  essentially  terminating  around  147  ns  attaining 
a  peak  value  at  145  ns.  Finally,  the  total  integrated  K-shell 
yield  is  shown  in  Fig.  8.  The  total  K-shell  yield  is  about  2.1 
kj  and  is  approximately  the  same  as  the  majority  of  GAMBLE  II 
results.  Maximum  K-shell  yield  for  GAMBLE  II  neon  experiments  is 
predicted  to  be  about  4  k j  .  This  optimized  yield  has  been 
achieved  experimentally.^ 

A  similar  set  of  numerical  simulations  were  carried  out  for 
the  case  of  the  PEOS.  The  initial  mass  per  unit  length  and  radii 
are  the  same  as  the  case  above.  However,  as  already  noted  the 
current  trace  is  different.  With  the  PEOS  the  current  waveform 
is  represented  by  the  profile  shown  in  Fig.  9.  The  peak  current 
is  about  0.88  MA,  sharper  current  risetime  and  a  broader  profile 
width.  The  current  also  remains  on  for  a  longer  period  of  time  - 
until  around  220  ns  -  a  full  85  ns  longer  than  the  no-switch 


case.  The  details  of  the  implosion  dynamics  are  shown  in  Figs. 
10-16.  Note  that  with  reduced  peak  current  it  takes  longer  to 
drive  the  plasma  in  and  results  in  a  somewhat  softer 
compression.  This  is  reflected  by  the  slower  implosion  velocity 
and  correspondi ng  cooler  temperatures  and  lower  densities  leading 
to  reduced  X-shell  yields.  These  results  indicate  that  peak 
current  and  its  maintainence  appear  to  be  the  only  parameters 
that  influence  significantly  the  yield,  assuming  some 
optimization  of  the  kinematics.  The  influence  of  the  current 
risetime  will  be  discussed  below.  This  revelation  is  not 

ii 

surprising  and  is  consistent  with  our  earlier  findings  and  is  in 
agreement  with  most  0-D  codes. 

We  will  now  discuss  the  results  of  a  series  of  numerical 
simulations  illustrating  the  systematic  behavior  of  the  kinematic 
variables  and  yield  as  a  function  of  mass  and  radius.  The  first 
series  of  simulations  varies  the  mass/length  while  maintaining 
the  initial  radii  fixed,  i.e.,  R Q  »  1.55  cm  and  =  0.95  cm  and 
M/ 2,  varies  between  15  to  35  pgm/cm  in  5  ugm/cm  steps.  Figure  17 
shows  the  maximum  velocity  as  a  function  K/2,  for  both  the  switch 
and  no-switch  cases.  For  a  fixed  peak  current,  the  higher  the 
M/2,  the  smaller  the  peak  velocity  of  implosion- pushing  on  the 
plasma  with  a  greater  magnetic  force  produces  a  greater  plasma 
acceleration.  These  results  suggest  the  following  scenario:  the 
higher  the  peak  current  the  higher  the  peak  temperature,  density, 
and  K-shell  yield,  provided,  of  course,  that  the  current  is  not 
so  high  that  it  burns  through  the  K-shell,  which  would  result  in 
a  reduced  K-shell  yield.  This  is  shown  in  Figs.  18-20  which 
compares  the  switch  and  no-switch  results.  The  higher  the  peak 
current  the  better  the  quality  of  the  implosion,  assuming 
optimized  kinematics,  again  provided  there  is  sufficient  mass  for 
good  K-shell  emission.  The  neon  K-shell  yield  presented  in  Fig. 
20  shows  that  below  about  22  ygm/cm  the  K-shell  yield  is  actually 
greater  for  the  lower  peak  current  (switch)  case  than  the  no¬ 
switch  case.  The  explanation  is  that  the  higher  current  no¬ 
switch  case  burns  through  much  of  the  K-shell  for  low  mass  loads 
and  hence,  has  fewer  K-shell  emitters  and  therefore  produces  a 


reduced  yield.  Above  22  ygm/cm  the  no-switc'n  yield  continues  to 
increase  as  the  available  number  of  K-shell  emitters  increases. 
Also  note  that  in  support  of  this  behavior  the  yields  of  each 
case  is  optimum  for  different  masses:  the  no-switch  case  at 
about  26  ugm/cm  and  the  PEOS  case  at  about  22  ugm/cm.  The 
percent  ratio  in  peak  current  is  virtually  the  ratio  of  yields; 
i.e.,  peak  currents  are  in  the  ratio  of  3/12.5  -  64$  while  the 
peak  K-shell  yields  are  in  the  ratio  of  2. 1/3.2  »  65$.  This  is 
not  meant  to  suggest  that  the  relationship  between  current  and 
yield  is  a  linear  one,  because  it  is  not,  it  is  only  indicative 
that  in  general,  for  optimum  kinematics,  the  higher  yields  are 
associated  with  higher  currents.  In  fact,  a  scaling  law  has  been 
derived  and  suggests  that  the  K-shell  yield  scales  approximately 
as  the  fourth  power  of  the  peak  current  (actually  1^“®). 

The  influence  of  current  risetime  on  the  yield  is  not  as 
obvious  from  these  calculations  as  it  is  from  the  experiments. 
The  experiments  indicate  a  sharper  risetime  current  profile,  and 
larger  yield  of  the  radiation  K-shell  pulse  when  the  PEOS  is 
used.  Also,  the  plasma  appears  to  exhibit  less  flaring  and  more 
uniformity  with  the  PEOS.  C al cul at i onal 1 y ,  we  cannot  comment  on 
the  plasma  structure  with  the  0-D  code  but  the  0-D  results 
predict  good  yields.  In  fact,  a  series  of  runs  were  made  with  a 
reduced  peak  current  in  the  no-switch  case  equal  to  the  peak 
current  in  the  PEOS  case.  The  yield  was  half  the  PEOS  yields 
supporting  the  notion  that  sharper  current  risetime's  influence 
the  yield.  It's  as  if  the  current  "hits"  the  load  so  quickly 
that  the  load  just  percolates  and  then  coasts  inward  producing  a 
more  uniform  implosion.  Obviously  this  is  an  area  that  requires 
further  investigation  if  we  are  to  understand  the  influence  of 
sharp  current  risetime  pulses  on  the  implosion  processes. 

The  final  set  of  numerical  simulations  involves  keeping  the 
mas3  per  unit  length  fixed  while  varying  the  initial  sheath 
thickness.  The  conditions  are  for  M / l »3 Ougm/ cm  and  a  fixed  outer 
radius  of  1.55  cm .  The  results  of  the  calculations  are  shown  in 
Fig3.  21-24  for  both  current  cases.  The  peak  velocities  decrease 
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slightly  with  increasing  Ar  which  just  reflects  the  behavior  of 
the  outer  boundary  vis  a  vis  the  inner  boundary:  i.e.,  since  the 
mass,  outer  radius,  and  current  are  fixed  the  only  variable 
parameter  is  the  inner  radius.  Depending  on  the  time  it  takes 
the  inner  radius  to  collapse  influences  the  competition  between 
magnetic  and  fluid  pressure  and  the  bounce  time  of  the  outer 
boundary.  In  Fig.  22  the  peak  densities  are  plotted  as  a 
function  of  Ar .  The  higher  densities  occur  for  lower  values  of 
Ar  and  result  from  more  particles  per  unit  volume  initially  (even 
though  the  total  number  of  particles  per  gram  is  the  same  in  all 
cases)  combined  with  a  better  implosion,  i.e.,  higher  velocity 
and  greater  compression.  For  larger  Ar  the  difference  between 
the  switch  and  no-switch  cases  is  not  readily  disoernable  on  the 
graph  but  the  no-switch  case  is  about  a  factor  of  two  greater  for 
the  higher  values  of  Ar .  The  peak  temperatures  as  a  function  of 
Ar  are  shown  in  Fig.  23.  At  first  glance  the  results  appear 
co  nt  r  adi  ct  or  y  in  the  sense  that  it  might  be  anticipated  that  the 
temperature  would  decrease  as  the  sheath  thickness  increased 
since  the  implosion  exhibited  slower  peak  velocities  and  lower 
peak  densities.  Actually,  the  peak  temperatures  increase  because 
there  are  fewer  emitters  which  reduces  the  major  source  of 
radiation  cooling  -  line  radiation  -  keeping  the  plasma  hotter. 
This  is  further  evidenced  in  Fig.  24  where  the  peak  radiative 
yields  are  shown  as  a  function  of  Ar ,  i.e.,  the  yield  (or 

alternatively,  the  radiation  cooling)  decreases  drastically  for 
large  Ar.  Since  the  yield  is  shown  only  for  the  K-shell,  a 

logical  conclusion  is  that  the  temperature  becomes  hot  enough  to 
erode  the  K-shell,  i.e.,  reduce  the  number  of  emitters,  and 
reduce  the  plasma  cooling  causing  the  plasma  to  remain  hot. 
Finally,  the  radiation  pulses  emitted  by  the  plasma 
experimentally  are  longer  than  those  calculated.  This  effect  is 
probably  due  to’  the  "zippering"  effect  of  the  implosion  and  is 
not  included  in  the  model.  Efforts  to  stagger  the  implosion  to 
simulate  this  effect  are  currently  underway. 


A  series  of  numerical  simulations  were  performed  to 
investigate  the  behavior  of  an  imploding  neon  gas  puff  for 
conditions  similar  to  experiments  conducted  on  the  NRL  GAMBLE  II 
pulsed  power  facility.  The  calculations  were  done  using  the 
SIMPLODE  code  which  is  best  described  as  a  0-D  Non-LTE  Dynamic 
Bennet  Pinch  model.  The  Non-LTE  model  was  based  on  a 

col 1 i s i onal -rad i a t i ve  equilibrium  description  in  the  optically 
thin  approximation.  Line  radiation  from  the  K-shell  only  was 
presented  although  total  radiation  was  accounted  for  in  the 
radiative  cooling  term.  The  results  of  the  numerical  simulations 
are  in  reasonably  good  agreement  with  the  experiments  -  at  least 
to  the  extent  of  the  limited  comparison.  As  more  experimental 
data  is  unfolded,  particularly,  emission  spectra,  a  more  thorough 
analysis  should  lead  to  a  better  understanding  of  how  to  control 
these  implosions.  What  we  have  learned  from  these  simulations  is 
that  impressive  K-shell  yields  can  be  obtained  with  neon  on  the 
GAMBLE  II  facility  and  that  peak  current  (and  possibly  current 
risetime),  assuming  optimized  kinematics,  is  the  single  most 
important  parameter  in  influencing  radiative  yields. 
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